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Abstract 

A non-LTE analysis of C i lines at 1.068-1.069 was carried out for selected 46 halo/disk stars covering 
a wide metallicity range (—3.7 ^ [Fe/H] ^ +0.3), based on the spectral data collected with IRCS+A0188 of 
the Subaru Telescope, in order to examine whether and how these strong neutral carbon lines of multiplet 
1 can be exploited for establishing stellar carbon abundances, especially for very metal-poor stars where 
CH molecular lines have been commonly used. These C i lines were confirmed to be clearly visible for all 
stars down to [Fc/H] ~ —3.7, from which C abundances could be successfully determined. The resulting 
[C/Fe] vs. [Fe/H] diagram revealed almost the same trend established from previous studies. When the 
results for individual stars are compared with the published data collected from various literature, while 
a reasonable agreement is seen as a whole, a tendency is observed that our abundances arc appreciably 
higher than those from CH lines especially for very metal-poor giants of low gravity. Since the abundances 
of these C i lines are subject to rather large non-LTE corrections (typically by several tenths dex) whose 
importance progressively grows as the metallicity is lowered, attention should be paid to how the coUisional 
rates (especially due to neutral hydrogen) are treated in non-LTE calculations. 

Key words: stars: abundances — stars: atmospheres — stars: late-type 
- stars: Population II 



1. Introduction 

The abundances of carbon (one of the most abundant 
elements next to hydrogen, helium, and oxygen) play a 
significant role in stellar spectroscopy. Above all, the be- 
havior of C in the very metal-poor regime ([Fe/H] ^ —2) 
is important, since various peculiarities (i.e., enhancement 
as well as deficiency) are observed, which have attracted 
interest of astrophysicists in the field of galactic nucle- 
osynthesis or stellar evolution. Thus a number of carbon 
abundance studies of metal-deficient stars have been pub- 
lished so far, as compiled in SAGA^ (Stellar Abundances 
for the Galactic Archeology) database (Suda et al. 2008, 
2011). Regarding the spectral lines to be invoked for spec- 
troscopically determining the abundance of carbon, sev- 
eral indicators are known to usable, such as atomic lines 
(C i) or molecular lines (CH, C2), each of which have dif- 
ferent visibilities depending on the type of stars, and have 
different characteristics; e.g., in terms of the sensitivity 
to the 3D effect or non-LTE effect (cf. Asplund 2005; 
subsection 3.3 therein). 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 
The large data tables are separately provided in the machine- 
readable form as electronic tables El and E2. 
Available on-line at (http://saga.sci.hokudai.ac.jp/). 



Focusing on very metal-deficient stars in the galactic 
halo (—4 ^ [Fe/H] ^ —2), we note that most studies on 
carbon abundances have been done based on CH molecu- 
lar lines in the blue region (G band) which are so strong 
as to be visible even at such an extremely low metallicity. 
Although a group of neutral carbon lines at 9061-9111 A 
(multiplet 3; 3s - 3p 

; Xiow ^ 7.5 eV) were some- 
times used for metal-poor dwarfs by several investigators 
(Tomkin et al. 1992; Akerman et al. 2004; Takeda & 
Honda 2005; Fabbian et al. 2006, 2009), they have not 
been the mainstream (presumably because they may suf- 
fer an appreciable non-LTE effect). 

Here, we should realize that other C i lines are available 
in near-IR region, which are namely the lines at ~ 1.07 /xm 
(multiplet 1; 3s - 3p ^D; 

Xiow '^7.5 cV) originating 
from the same lower term (3s ^P°) as multiplet 3 lines. 
Actually, they are even more suitable than C i 0.91 /im 
lines, because of their strengths (the transition probability 
of the strongest C i 10691 line in multiplet 1 is larger by 
0.2 dex than the strongest C I 9094.83 line in multiplet 
3; cf. table 2) as well as their location in the spectral 
region almost free from telluric lines (these are actually a 
nuisance problem for C i 0.91 fj,m lines). However, to our 
knowledge, these lines have barely been used for carbon 
abundance determinations, except for the Sun and Vega 
(e,g., Stiirenburg & Holweger 1990; Takeda 1992, 1994). 
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Why not exploit these lines to study carbon abundances 
of mctal-poor stars by making use of these merit? 

Timely, we recently obtained near-IR (zJ-band; 1.04- 
1.19 ^m) high-dispersion spectra of 46 disk/halo stars 
(dwarfs as well giants) in a wide metallicity range 
(-3.7 < [Fe/H] < +0.3), by using the Subaru Telescope 
with IRCS+A0188, and carried out sulfur abundance 
determinations based on the S i triplet lines at 10455- 
10459 A (Takeda & Takada-Hidai 2011a, 2012; hereinafter 
referred to as Paper I and Paper II, respectively). We 
thus decided to study the strong C i lines of multiplet 1 
at 1.07 fiia in these spectra, in order see whether the 
C abundances of very metal-poor stars can be success- 
fully determined by these lines and how they are compared 
with the published results derived from other lines (such 
as those of CH) . This is the purpose of this investigation. 

The remainder of this article is organized as follows. 
After describing the observational data in section 2, we 
explain the details of abundance determinations (stellar 
parameters, treatment of non-LTE effect, method of anal- 
ysis, uncertainties, etc.) in section 3. Section 4 is devoted 
to the discussion section, where the results are examined, 
compared with the published work, and the merits and 
problems of these near-IR C i lines are discussed. The 
summary is given in section 5. In the Appendix, the grid 
of theoretical non-LTE corrections for the C i 9061-9111 
lines of multiplet 3 are re-presented, since a mis-treatment 
was found in our previously published results. 

2. Observational Data 

The observational data used for this study are the high- 
dispersion (i? ^ 20000) zJ-band (1.04-1.19 ^m) spectra 
obtained by using IRCS+A0188 of the Subaru Telescope 
in two observing periods of 2009 July 29 and 30 (UT) as 
weh as 2011 August 17 and 18 (UT). In the former 2009 
July run, 33 comparatively bright stars of wide varieties 
(from [Fe/H] ~ —3.7 to ^ -1-0.3, dwarfs as well as giants) 
were observed,^ while selected 13 very metal-poor stars 
([Fe/H] < —2) and Vesta (substitute for the Sun) were 
targeted in the latter 2011 August run. This makes 47 
(= 33-1-13-1-1) spectra as a total, though the net number is 
46 because G 64-37 was repeatedly observed. See section 
2 of Papers I and II for more details of the observations 
and the data reduction. The list of these 47 objects is 
given in table 1, where the targets in 2009 and those in 
2011 are separately presented. 

The 2009 July spectra for representative 4 stars of differ- 
ent metallicities along with Altair in the 10650-10750 A 
region are displayed in figure 1. Note that most of the 
conspicuous lines are those of neutral carbon and sili- 
con, and several strong C i hnes (e.g., at 10683, 10691, 
10730 A) are clearly visible also for the most metal-poor 
star (BD -1-44 493). Besides, we can recognize from this 
figure (see the broad-line spectrum of Altair) that this re- 
gion is practically free from any telluric lines, in contrast 

^ In addition, Altair (A7 V, rapid rotator with Vc sin i ~ 
200 km s~^) was also observed to be used as a reference of 
telluric lines. 
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Fig. 1. Spectra of representative 4 stars with different metal- 
licities in the 10650-10750 A region comprising C I and Si I 
lines, where the wavelength scale of stellar lines is adjusted to 
the laboratory frame. The spectrum of Altair is also shown 
to demonstrate the negligible effect of telluric lines. 

to the C I 9061-9111 region influenced by H2O lines of 
earth's atmosphere origin. 

3. Abundance Determination 

3.1. Non-LTE Calculations and Model Atmospheres 

Regarding the atmospheric parameters (Tes, \ogg, Vt, 
and [Fe/H]) of the program stars necessary for construct- 
ing model atmospheres and determining abundances, we 
used exactly the same values as adopted in Papers I and 
II, in which various published studies were consulted (cf. 
caption in table 1 of Papers I and II for the individual ref- 
erence sources). These parameters are presented in table 
1. 

As to the computation of non-LTE departure coeffi- 
cients, we closely followed our previous work (see Takeda 
1992, 1994; Takeda & Honda 2005 for the details). 
Practically, our non-LTE calculations were carried out on 
a large grid of 150 (=6x5x5) models corresponding to 
combinations of six Teff (4500, 5000, 5500, 6000, 6500, 
7000 K), five \ogg (1, 2, 3, 4, 5), and five [Fe/H] (-4, -3, 
—2, —1, 0), where [C/Fe] = (carbon abundance) and vt 
~ 2 km s~^ (microturbulence) were assumed. 

To keep consistency with our previous studies, the coU- 
sional rates (Co for electron collisions and Ch for neu- 
tral hydrogen collisions) were computed by following the 
recipe described in subsubsection 3.1.3 of Takeda (1991). 
Namely, regarding electron collisions, we invoked Van 
Regemorter's (1962) formula for permitted transitions 

{C'^^) and Auer and Mihalas's (1973) equation (13) for 
forbidden transitions (Cc ). Meanwhile, as for hydrogen 
collision rates, for permitted transition was evaluated 
by invoking Steenbock and Holweger's (1984) formula, 
which is based on Drawin's (1968, 1969) classical cross 
section, and for forbidden transitions was approxi- 

mated by simply scaling Ce according to the collision 
frequency (determined by the relative speed and the den- 
sity of colliding particles) on the assumption of the same 
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cross section. Since such computed Ch(= +C'h^) has 
only ordcr-of-magnitude accuracy at best and tends to 
be often overestimated, we also tested an alternative case 
of reduced Ch by introducing a reduction factor k as was 
done in Takeda and Honda (2005; cf. appendix 1 therein): 

Ctotal=Ce + fcC7H, (1) 

though we still assume fc = 1 as the standard case to be 
adopted. Such a scaling of Ch with k is applied also to 
bound-free collisional ionization rates, for which we in- 
voked Jefferies (1968) for electron collisions and Steenbock 
and Holweger (1984) for neutral- hydrogen collisions. 

We then interpolated Kurucz's (1993) grid of ATLAS9 
model atmospheres as well as the grid of the non-LTE de- 
parture coefficients in terms of Toff, \ogg, and [Fe/H] to 
generate the atmospheric model and the departure coeffi- 
cient data for each star. 

Besides, for the reader's convenience, extensive grids of 
non-LTE as well as LTE equivalent widths {EW^ and 
EW^) and non-LTE abundance corrections (A) for five 
multiplet 1 lines (C i 10683, 10685, 10691, 10707, 10729) 
were further computed for three [C/Fe] values (—0.3, 0.0, 
and -1-0.3) and three Vt values (1, 2, and 3 km s~^), based 
on the non-LTE departure coefficients computed for 150 
model atmospheres with two different treatments of neu- 
tral hydrogen collisions (k — 1 and k ~ 0.1), which are 
presented in electronic table El. 

3.2. Synthetic Spectrum Fitting 

Inspecting the spectrum feature around ~ 1.07 fj.m (cf. 
figure 1), we decided to concentrate on the 10680-10697 A 
region which comprises three C i lines (in which C i 10691 
is the strongest one of multiplet 1) and two Si i lines. 
Then, following the similar way as in Papers I and II, we 
carried out non-LTE spectrum-synthesis analyses by ap- 
plying Takeda's (1995) automatic fitting procedure to this 
17 A-interval region while regarding the non-LTE carbon 
abundance (A^), LTE silicon abundance (Agj), the macro- 
broadening parameter (wm; e- folding half- width of the 
Gaussian macrobroadening function; c>c exp[— (u/tjm)^]), 
and the wavelength shift (AA) as adjustable free param- 
eters to be varied until a convergence is accomplished. 
The adopted atomic data of the relevant C i and Si i lines 
are presented in table 2. This fitting procedure turned 
out quite successful and could be established for all 
47 cases, though we had to fix the silicon abundance for 
two extremely metal-poor stars (BD +44 493 and G 64-37 
[2009 July data]) where Si i lines are hardly visible. How 
the theoretical spectrum for the converged solutions fits 
well with the observed spectrum is displayed in figure 2 
(2009 July data) as well as figure 3 (2011 August data). 
(Note that these two figures are arranged in analogy with 
figure 2 of Paper I and figure 1 of Paper II, respectively.) 
The resulting non- LTE C abundances {A^) are given m 
table 1, where the values of [C/Fe] {=A^-A^q) - [Fe/H]; 
with Aq =8.54 from Vesta) are also presented. More de- 
tailed results are found in electronic table E2, where the 
results of and dm are also given. 



3.3. Abundance- Related Quantities 

While the non-LTE synthetic spectrum fitting directly 
yields the final abundance solution, this approach is not 
necessarily suitable when one wants to evaluate the ex- 
tent of non-LTE corrections or to study the abundance 
sensitivity to changing the atmospheric parameters (i.e., 
it is tedious to repeat the fitting process again and again 
for different assumptions or different atmospheric pa- 
rameters). Therefore, with the help of Kurucz's (1993) 
WIDTH9 program (which had been considerably mod- 
ified in various respects; e.g., inclusion of non-LTE ef- 
fects, etc.), we computed the equivalent widths for each 
of the three C i lines (£^M^io683, EWwesb, and EWwegi) 
"inversely" from the abundance solution (resulting from 
non-LTE spectrum synthesis) along with the adopted at- 
mospheric model/parameters, since they are much eas- 
ier to handle. Based on such evaluated EW values, the 
LTE abundances for each of the lines (^loess i ^loess ' ^^'^ 
^10691 ) were freshly computed, from which the non-LTE 
corrections (Aioess, Aioess, and Aioggi) were derived such 

as Aio683 = ^10683 " ^C' ^tC. 

Besides, in order to examine the effect of changing 
Ch (neutral hydrogen collision rates), we also computed 
Aio69i {k — 0.1) from EWioegi by using the non-LTE de- 
parture coefficients calculated for the fc = 0.1 case (i.e., 
Ch is reduced by 1/10 from the adopted standard case of 
fc = 1) and evaluated the difference 

(5Aio69i = Aio69i(A: = 0.1) - Aio69i(fc = 1). (2) 

We then estimated the uncertainties in A^ by repeat- 
ing the analysis on EWioe9i {EW for the strongest com- 
ponent) while perturbing the standard values of atmo- 
spheric parameters interchangeably by ±100 K in Toff, 
±0.2 dex in logg, and ±0.3 km s~^ in vt (which we re- 
garded as typical uncertainties of the atmospheric param- 
eters according to the original references; cf. subsection 
4.2 in Paper I). Let us call these six kinds of abundance 
variations as St+, 5t-, Sg+, Sg-, Sv+, and Sy-, respec- 
tively. We then computed the root-sum-square of three 
quantities Sxgv = (i^t + ±'^^)^^^ as the abundance un- 
certainty (due to combined errors in Tcs, logg, and Vt), 
where St, Sg, and Sy are defined as St = (|i5t+| ± |(5t-|)/2, 

Sg = {\Sg+\ + \Sg-\)/2, Rlid Sy = {\Sy+\ + \Sy-\)/2, rCSpCC- 

tively. 

The resulting [C/Fe], Aio69i, <5Aio69i, and EWioegi for 
each star are plotted against Toff in figures 4 a-d, where 
logg > 3.0 stars (dwarfs) and logg < 3.0 stars (giants) are 
distinguished by filled and open symbols, respectively, and 
the error bar attached to [C/Fe] represents ±5Tgy While 
only A^, EWioQQi, Aio69i, and [C/Fe] are given in table 
1, all the relevant data (including EW and A for all three 
lines and each of the S values) are presented in electronic 
table E2. 

It can be seen from figure 4a that the typical extent 
of STgy is ^ 0.1-0.2 dex and not very significant. We 
should bear in mind, however, that the non-LTE correc- 
tions are so large and that some uncertainties in them are 
more or less inevitable. Apart from the ambiguities in 
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Fig. 4. Carbon abundances (along with the related quan- 
tities) plotted against [Fe/H]: (a) [C/Fe] (C-to-Fe logarith- 
mic abundance ratio corresponding to non-LTE carbon abun- 
dance; where attached error bars represent the ambiguities 
due to uncertainties in the atmospheric parameters {5Tgv\ 
cf. subsection 3.3) (b) Aioegi (non-LTE correction for the 
C I 10691 line derived with the standard classical treatment 
of neutral-hydrogen collisions with k = 1). (c) 5Aio69i (rela- 
tive difference of non-LTE correction caused by reducing the 
neutral-hydrogen collisions by 1/10; cf. equation (2)). (d) 
^'M^10691 (equivalent width for the C I 10691 line). Dwarfs 
(logg > 3) and giants (logg < 3) are discriminated by filled 
(blue) and open (green) symbols, respectively. The results 
for 33 stars based on the 2009 July data and those for 12 
stars based on the 2011 August data are shown by circles and 
squares, respectively. The large red filled circle denotes the 
Sun (Vesta). 

the treatment of coUisional rates due to neutral hydrogen 
collisions (choice of k) as discussed later in subsections 
4.1 and 4.3, we should recall that our non-LTE calcula- 
tions were done on the condition of [C/Fe] = 0. Actually, 
this assumption is reasonable, since most of the [C/Fe] 
values derived for our targets scatter around with a 
dispersion of ~ ±0.5 dex (figure 4a). However, [C/Fe] val- 
ues of two stars considerably depart from zero; i.e., -I-0.83 
(BD -t-44 493) and -0.94 (HD 13979), for which the mis- 
match of [C/Fe] may cause some errors. As a test, we 
evaluated the abundance change when we used departure 
coefficients computed with [C/Fe] = +1.0 (BD +44 493) 
and [C/Fc] = -1.0 (HD 13979) instead of [C/Fe] = 0. 
These calculations were done on the selected two models 
(t55g40m40 and t50g20m20 included in the model grid of 
electronic table El) having parameters closest to those 
of these two stars. We then found that the extent of 



non-LTE overpopulation for the relevant high-excitation 
term is enhanced/reduced by decreasing/increasing the C 
abundance,"^ and thus the Aq (derived for the same EW) 
is increased by 0.25 dex ([C/Fe] = +1.0; BD +44 493) 
and reduced by 0.10 dex ([C/Fe] = -1.0; HD 13979) 
compared to the results for the standard [C/Fe] = 0.0 
case. Accordingly, the error caused by applying the de- 
parture coefficients computed with [C/Fe] = to such 
exceptional cases of considerably peculiar C abundance 
would be 0.1-0.2 dex. 

4. Discussion 

Line Strength and Non-LTE Effect 

According to figures 4b-d and table 1 , the following fea- 
tures are recognized regarding the behavior of the strength 
and the non-LTE correction of the C i 10691 line (the 
strongest among the multiplet 1 lines at A 1.07 ^m). 

— The C I 10691 line is clearly visible even for extremely 
metal-poor stars ([Fe/H] < —3) with EW on the order of 

10 mA, irrespective of dwarfs or giants, which suggests 
the usability of this line to explore the carbon synthesis 
history of the Galaxy. 

— The non-LTE correction is always negative (corre- 
sponding to non-LTE line-strengthening). While the typi- 
cal extent of [Aioegil for the standard fc = 1 case is several 
tenths dex (^ 0.5 dex; cf. figure 4b), it is further raised 
by almost a similar amount when k is reduced from 1 to 
0.1 (figure 4d) which makes [Aio69i(fc ~ 0.1) | even up to 
^ 1 dex. This implies that the non-LTE correction is sen- 
sitive to how the neutral hydrogen collision is treated.^ 

— The extent of |Aio69i| systematically increases with a 
decrease in [Fe/H], despite that the line strength is weak- 
ened, which means that the non-LTE correction becomes 
progressively more important as we go toward the ex- 
tremely metal-poor regime. 

In order to help our understanding of these features 
mentioned above, we display in figure 5 how EW^ (the- 
oretical non-LTE equivalent width computed for the C i 



This trend could be interpreted as due to the efficiency change 
in the UV pumping from ground or low-excitation levels to up- 
per levels of high excitation, which may possibly depend upon 
the C abundance in the very metal-poor regime. That is, as 
long as such UV transitions are optically thick and in radiative 
detailed balance, such process does not occur. Meanwhile, as 
the C abundance is sufficiently lowered and some of these UV 
transitions eventually become optically thin, high-excitation C I 
levels may get overpupulated via upward pumping by absorbing 
hot UV continuum radiation. 

This does not agree with the conclusion of Fabbian et al. (2006, 
2009), who found that the non-LTE corrections for C I lines 
are insensitive to a choice of the reduction factor (Sjj in their 
notation); even changing Sh = 1 to 5h = makes only an in- 
significant change (~ 0.1 dex) in their non-LTE corrections for 
C I 9061-9111 lines of multiplet 3, which are as large as those of 
C I 10683-10691 lines of multiplet 1 (cf. figure 8b). Considering 
the more or less reasonable consistency between the A (multiplet 
3) values of Takeda and Honda (2005) (see also the Appendix 
of this paper) and those of Fabbian et al. (2006; cf. their figure 
8), this discrepancy is hard to understand. Although the cause 
is not clear, something may be different between their and our 
calculation procedures. 
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Fig. 5. Theoretically computed non-LTE logarithmic equiv- 
alent widths (log EW^; lower panels) and the difference 
between non-LTE and LTE logarithmic equivalent widths 
(log EW^ — log EW^; upper panels) plotted against T^g, 
which were calculated for the C I 10691.24 line with the stan- 
dard treatment of H I collision cross sections (fe = 1), vt = 
2 km s~^, and [C/Fc] = 0. The left, middle, and right panels 
correspond to [Fe/H] = 0, —2, and —4, respectively. Open 
circles, filled circles, open triangles, filled triangles, and open 
squares denote results for logg = 1, 2, 3, 4, and 5, respectively. 

10691 line) and EW^/EW^ (non-LTE to LTE equivalent- 
width ratio) depend on Tcff, logg, and [Fe/H], based on 
the data in electronic table El (cf. subsection 3.1) corre- 
sponding to the standard k ^ 1 case. We can learn from 
this figure significant characteristics concerning the line- 
strength behaviors: 

— While it is natural that the line is weakened with a 
decrease in [Fe/H] (=[C/H]), its strength grows with an 
increase in Toff as well as with a decrease in logg. 

— The non-LTE to LTE strength ratio is always larger 
than unity (i.e., non-LTE line strengthening) as we al- 
ready recognized from the negative sign of A, which in- 
creases with a lowering of the mctallicity; it amounts up to 
^5 at [Fe/H] = —4, which explains the [Fe/H]-dependence 
of |A| mentioned above. 

— We may interpret these trends in terms of the ioniza- 
tion degree of carbon (partially ionized in the atmospheres 
of FGK stars), which tends to enhance with an increase in 
Toff as well as with a decrease in logg (lowered density) and 
[Fe/H] (lowered electron density, enhanced UV overioniz- 
ing radiation due to the increased transparency), while 
the change in the H~ continuum opacity is also involved 
especially in the logg-effect. That is, the population of 
high-excitation C i lines is closely coupled with the popu- 
lation of C II parent term (i.e., the effect of n^^^'f /n^^'>^; cf. 
section 3 in Takeda 1992), which makes the population of 
high-excitation C i level larger as the ionization fraction 
(C ii/C i) increases. Therefore, the degree of ionization 
plays an important role in understanding the strengths 
and the non-LTE effect of high-excited C I lines such as 
those of multiplet 1 or 3. 

— Besides, under the considerably metal-poor condition 



of very low C abundance, some of the strong C i transi- 
tions in UV would become optically thin, which may act to 
enhance non-LTE overpopulation of high-excitation levels 
via the pumping effect (cf. footnote 3). 

— Thus, the mechanism of non-LTE line intensification is 
not necessarily the same for the strong-line case of near- 
solar metallicity and the weak-line case of very low metal- 
licity. In the former situation, the line is strengthened sig- 
nificantly by the dilution of the line source function due 
to photon escape (such as the case of a simple two-level 
atom), while the increase of line-opacity due to overpop- 
ulation in the lower level of this transition is essential for 
the latter case. 

4-2. Comparison with Literature Data 

We now discuss the carbon abundances we have derived 
in comparison with previous studies taken from the liter- 
ature. 

According to figure 4a, our [C/Fe] vs. [Fe/H] relation 
constructed from 47 objects is characterized by (i) a grad- 
ual increase of [C/Fe] at a slightly supersolar level with a 
decrease of [Fe/H] from ~ to — 1, (u) a subsolar [C/Fe] 
around ~ —0.5 (mostly giants) at —3 ^ [Fe/H] ^ —2, and 
(iii) upturn of [C/Fe] toward appearance of supersolar 
[C/Fe] at the extremely metal-poor regime (— 4< [Fe/H] ^ 
—3). Actually, these features arc in good agreement with 
what has been established so far; e.g., figure 6 of Norris, 
Ryan, and Beers (1997) is quite similar to our figure 4a. 
Accordingly, we may state that our analysis using C i 
10683/10685/10691 lines has yielded [C/Fe] results con- 
sistent with other previous work (where CH lines were 
mainly used in the very metal-poor regime) at least in the 
qualitative sense. 

We then turn to the comparison of our [C/H] (= Aq — 
Aqq ) for individual stars with the literature data summa- 
rized in table 3, which we collected mainly by invoking the 
SAGA database (cf. section 1) and by adding data from 
several papers by ourselves (see the caption of table 3). 
We could thus get the published data for 37 stars (^ 80% 
of the total objects). These literature [C/H] values are 
plotted by crosses (red • • • molecular lines mostly of CH, 
blue • • • C I lines) in figure 6, where our results are desig- 
nated by circles (open ■ ■ ■ giants, filled • • • dwarfs). Several 
notable characteristics can be read from this figure: 

— For disk stars with metallicity of — 1 < [Fe/H] where 
C I lines are mostly used, the agreement is generally good, 
without any systematic discrepancy. 

— In contrast, at the very low mctallicity range 
(-2.7< [Fc/H] < -1), we note that literature [C/H] values 
derived from CH lines (red crosses) tend to be appreciably 
lower (typically by several tenths dex) than our results, es- 
pecially for low-gravity giants (open circles).^ Meanwhile, 



Actually, such a CH vs. C I discrepancy is apparently observed 
in the results of Tomkin et al. (1992), which is a rare study 
because two kinds of C abundances from C I lines and CH 
molecular band are given for a star. According to their table 
4, the average difference for 28 stars (in which both abundances 
are available) is (y4(CH) — A{C I)) = —0.41 dex, which is nearly 
on the same order as we see in figure 6. We should note, how- 
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Fig. 6. Comparison of [C/H] {= Ac — Aqq) values derived 
in this study (circles, where open and filled symbols denote 
giants and dwarfs, respectively) with those taken from various 
published studies (crosses, where those derived from molecu- 
lar and atomic lines are discriminated in red and blue colors, 
respectively), based on the data summarized in table 3. Stars 
are arranged in the ascending order of [Fo/H] from top to 
bottom. 

those derived from C i lines (mostly 9061-9111 lines of 
multiplet 3) expressed by blue crosses do not show such a 
systematic tendency. 

— However, if we confine to extremely metal-poor stars 
at [Fe/H] < —3, such a trend can not be clearly seen; e.g., 
CH-bascd literature [C/H]'s are even higher than our re- 
suhs in HD 126587 and BD +AA 493. 

— Roughly summarizing, we may state that our C- 
abundances determined from C i 10683-10691 lines are 
almost consistent with the literature results if they are 
based on C i lines, while CH-based abundances published 
so far tend to be lower than ours (especially for giants). 

4-. 3. Merits and Problems of C I 1.07fim Lines 

Returning to the motivation of this study described in 
section 1, we again question to ourselves: Are C i 1.07 jim 
lines of multiplet 1 practically useful for stellar carbon 
abundance determinations, especially for very metal-poor 

ever, that the real discordance could be much more serious (even 
amounting up to ^ 1 dex), since A(CH) may further be signifi- 
cantly lowered (e.g., by a few tenths dex to ^ 1 dex; especially 
at the very metal-poor condition) when new 3D model atmo- 
spheres are used (e.g., Asplund 2005; Collet 2007) because of 
the reduced temperature at the upper layers. 



stars? Our answer is as follows: "Yes, definitely. We 
believe that these near-IR lines would promisingly be ap- 
plied to investigate C abundances for a variety of stars. 
However, several things still remain to be done before we 
can assure a sufficient precision in abundance determina- 
tions." 

We summarize below the favorable characteristic of 
these C i lines in comparison to other indicators. 

— These C l lines have sufficiently large strengths and 
would be visible down to extremely low-metallicity re- 
gion. According to figure 5, if we could measure a very 
weak line of EW ^10 mA based on a high dispersion 
(e.g., R ~ 30000-50000) and high S/N (e.g., - 300-500) 
spectrum, Aq of a [Fe/H] ~ —A star would be deter- 
minable even for the case of [C/Fe] ~ (of course, de- 
tection becomes much easier for carbon-enhanced stars of 
[C/Fe] > 0). 

— This merit almost equally applies to stars of any type, 
since whether a star is a giant or a dwarf does not sig- 
nificantly affect the line strength. That is, since dwarfs 
with higher logg tend to have higher Tee (e.g., turn-off 
stars with Tgg ~ 6000 K and \ogg ~ 4) and giants with 
lower \ogg tend to have lower Teff (e.g., evolved red giants 
with Teff ^ 4000 K and log, 9 1) as shown in figure la of 
Paper I, these two opposite effect of Toff a-nd logg tend to 
cancel with each other (cf. figure 5), which makes the line 
strength rather insensitive to Toff as well as logg. 

— Carbon abundances from CH molecular lines are very 
sensitive to (cf. table 5 of Tomkin et al. 1992), as well 
as to the temperature structure in the upper layer and to 
the 3D effect (e.g., Asplund 2005; Behara et al. 2008). 
In this respect, C i lines arc much less sensitive to such 
effects. 

— These C i 10683-10691 lines of multiplet 1 are more 
advantageous than the C i 9061-9111 lines of multiplet 3, 
because they are stronger and situate in a spectral region 
almost free from telluric lines. 

Despite these merits, however, we have to keep in mind 
that problems still remain in deriving the C abundances 
from these C i lines, especially when a sufficient precision 
is pursued. Above all, given that these lines are subject to 
rather a large non-LTE effect, attention should be paid to 
use as reliable non-LTE abundance corrections as possible. 
Here, the difficult problem is the treatment of collisions 
with neutral hydrogen atoms, or how to choose the value 
of k in equation (1). Since the non-LTE corrections ap- 
preciably depend upon this parameter (as we can see from 
figure 4c), its adequate choice is mandatory. 

Yet, we consider that fc = 1 adopted in this investiga- 
tion (as in our past studies) is reasonable for the following 
reasons: 

— If we use k = 0.1, for example, we would obtain 
~ ^-"^0 (A = —0.28) as the solar carbon abundance 

(instead of 8.54 for k = 1). Although this value is almost 
the same as the most recent solar C abundance of 8.43 
claimed by Asplund et al. (2009) using the state-of-the- 
art 3D model, our result based on the classical ID model 
should not be directly compared to it. Considering that 
Takeda and Honda (2005) derived the solar C abundances 
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of 8.49, 8.76, and 8.70 (based on the same solar model as 
used in this study) from three non-LTE insensitive lines 
([C i] 8727, C I 5052, and C i 5380), we feel reluctant to 
further lower the value of 8.54. 

— Admittedly, the problematic discrepancy between our 
[C/H] results and those from CH molecular lines seen in 
metal-poor giants (cf. subsection 4.2) may significantly 
be mitigated by reducing k. If we do so, however, our 
C l-based abundances for dwarfs would be overcorrected 
(actually, non-LTE corrections for dwarfs are more sensi- 
tive to k than those of giants; cf. figure 4d) to become dis- 
cordant with the literature results, since the consistency 
between CH and C i results are not so bad in metal-poor 
dwarfs unlike the case for giants (cf. figure 6). We rather 
suspect that CH-based abundances of giants may have 
been incorrect. Given that modeling of giants atmosphere 
is generally less reliable than the case of dwarfs, previ- 
ously published C abundances of metal-poor giants based 
on CH lines may have been appreciably underestimated 
due to some imperfectness of the atmospheric model, re- 
flecting the enormous sensitivity of molecular lines to er- 
rors in the temperature structure. For example, since the 
ubiquitous existence of chromospheric activities even in 
very metal-poor old stars is evidenced by the detection 
of He I 10830 hue (cf. Takeda & Takada-Hidai 2011b), 
the temperature of molecule-forming upper layers may be 
higher than that predicted by classical model atmospheres 
or UV radiation coming from the chromosphere may pro- 
mote non-thermal dissociation of molecules in the pho- 
tosphere, which would lead to underestimation of carbon 
abundances derived from CH molecular band (if based on 
the conventional method of analysis) . 

In any event, it is very important to determine the value 
of k in an empirical manner. One way for doing this may 
be to require the abundance consistency between differ- 
ent lines, such as done by Takeda and Honda (2005, cf. 
appendix 1 therein) for oxygen lines. Actually we chal- 
lenged it by comparing the abundances from lines of mul- 
tiplet 1 and multiplet 3, as described in the Appendix. 
Unfortunately, this trial turned out unsuccessful, because 
the fc-sensitivity of these two line groups were essentially 
the same. As an alternative approach, it would be promis- 
ing to compare the center-to-limb variation of C i 10683- 
10691 lines on the solar disk with the simulations done for 
various k values, as recently tried by Pereira, Kiselman, 
and Asplund (2009a) as well as Pereira, Asplund, and 
Kiselman (2009b) for the oxygen lines. We would be able 
to use these near-IR C i lines for precise C abundance de- 
terminations with confidence, when the calibration of this 
coUisional parameter has been established. 

5. Conclusion 

Motivated by the fact that C i lines at ~ 1.07 fim (multi- 
plet 1) have barely been exploited for stellar C abundance 
studies so far, despite their sufficient strengths and be- 
ing in a favorable spectral region free from telluric lines, 
we carried out a non-LTE spectrum-fitting analysis of C i 
10683-19691 lines for selected 46 halo/disk stars in a wide 



7 

metallicity rang e (- -3.7 [Fe/H] < +0.3) based on the zJ- 
band spectral data collected with IRCS-I-A0188 of the 
Subaru Telescope, with an aim of examining whether the 
C abundances can be successfully determined by these 
lines especially for very metal-poor stars where lines of 
CH molecules have been mainly invoked. 

These C i lines were confirmed to be clearly visible, from 
which C abundances could be successfully determined for 
aU stars down to [Fe/H] -3.7. The resulting [C/Fe] 
vs. [Fe/H] diagram revealed the well-known trend estab- 
lished by previous studies; such as a gradual increase of 
[C/Fc] with a decrease of [Fe/H] from to ~ — 1, sub- 
solar [C/Fe] of 0.5 at -3 < [Fe/H] < -2, and up- 
turn of [C /Fe] toward appearance of supersolar [C /Fe] at 
—4 ^ [Fe/H] ^ —3. We may thus regard that our analysis 
is consistent with previous studies at least qualitatively. 

We also compared the resulting [C/H] for individual 
stars with the published data collected from various liter- 
ature (data for 37 stars were available). While the agree- 
ment is mostly good for disk stars (— 1 ^ [Fe/H]) where C i 
lines are generally used, we found that the literature [C /H] 
values of very metal-poor stars (—2.7^ [Fe/H] ^ —1) de- 
rived from CH lines tend to be appreciably lower than our 
results especially for low-gravity giants, though those de- 
rived from C I 9061-9111 lines of multiplet 3 do not show 
such a systematic tendency. The reason for this disagree- 
ment is yet to be investigated. 

Since the abundances of these neutral C lines are sub- 
ject to an appreciably large non-LTE effect (typically by 
several tenths dex) whose importance progressively grows 
as the metallicity is lowered, attention should be paid to 
the reliability of non-LTE corrections to be applied. Here, 
the difficult problem is how to treat collisions with neutral 
hydrogen atoms, for which classical rates are often multi- 
plied by a reduction factor k{< 1). Although we consider 
that the choice of fc = 1 (standard classical treatment) is 
still reasonable, its choice has a significant influence on the 
extent of non-LTE effect (e.g., if we use fc = 0.1, non-LTE 
corrections would become twice as large as the k = l case) . 
Thus, establishing fc in some empirical manner (e.g., by 
using the center-to-disk variation on the solar disk) is ur- 
gently awaited. 

One of the authors (M. T.-H.) is grateful for a flnancial 
support from a grant-in-aid for scientific research (C, No. 
22540255) from the Japan Society for the Promotion of 
Science. 

This research has made extensive use of the SAGA 
database system, maintained by Takuma Suda, Yutaka 
Katsuta, and Shimako Yamada. 

Appendix. Non-LTE Corrections of C I 9061- 
9111 Lines Revisited 

Takeda and Honda (2005; cf. appendix 2 therein) com- 
puted an extensive grid of non-LTE equivalent widths and 
abundance corrections for five C i lines of multiplet 3 
(3s 3p° - 3p 3p) at 9061.44, 9078.29, 9088.51, 9094.83, and 
9111.81 A. They also carried out a non-LTE reanalysis of 
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published equivalent- width data of these C i 9061-9111 
lines along with those of C i 7111-7119 lines (multiplet 
26, 25.02). 

Having reexamined those previous calculations, we re- 
alized that there was a mis-assignment of departure co- 
efhcicnts in computing the non-LTE equivalent widths of 
C I 9061-9111 lines. More precisely, we had erroneously 
used the departure coefficients corresponding to multiplet 
1 (3s '^P° - 3p '^D) for these multiplet 3 lines. Accordingly, 
we decided to redo the calculations (this time, computa- 
tions were performed not only for fc = 1 as assumed before, 
but also for A: = 0.1), and the resulting grids of new non- 
LTE corrections (A) and equivalent-widths (EW) are pre- 
sented in electronic table El along with those of C 1 10685- 
10729 hues. 

Comparing the new results with the old ones, we re- 
alized that the differences are insignificant; actually, the 
change of A is ^0.1 dex at most, as shown in figure 7. 
This is due to the fortunate fact that multiplets 1 and 3 
have the same lower term (3s '^P°), while the upper terms 
(3p "^D and 3p '^P) have similar departure coefficients be- 
cause they lie closely to each other. Accordingly, we do 
not see any need to apply significant changes to what is de- 
scribed in appendix 2 of Takeda and Honda (2005), which 
still remain essentially valid in the practical sense. 

In connection with the newly recomputed non-LTE cor- 
rections for multiplet 3 lines, wc tried to empirically de- 
termine the value of k defined in equation (1) by re- 
quiring the consistency between the carbon abundances 
derived from C i 10685-10691 lines (obtained from this 
study) and those from C i 9061-9111 lines (based on re- 
analysis of published equivalent widths). By consulting 
the work of Akerman et al. (2004), EWs of 9061-9111 
lines for 8 stars were available, from which we derived 
the mean abundance of each star for three choices of fc: 
k = 0.1, fc = 1, and k = oo (i.e., LTE). Comparison of such 
computed ^c(9061-9111:RMT 3) with our Ac(10685- 
10691 :RMT 1) is presented in figure 8a. Unfortunately, we 
can recognize from this figure that the fc-sensitivity of both 
multiplet lines is nearly the same with each other, which 
means that empirically establishing fc by the requirement 
of ^c(R.MT 1) = ^c(R.MT 3) is hardly practicable (cf. 
figure 8b). We would have to find some other different 
carbon lines to be compared, in order to accomplish this 
aim. 
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Fig. 7. (a) Comparison of the old non-LTE corrections 
(A°''^) for the 9094.83 line (the strongest one of the five 
9061-9111 lines of RMT 3) computed by Takeda and Honda 
(2005; cf. appendix 2 therein) for an extensive grid of model 
atmospheres, for which we found that departure coefficients 
for the upper term had been erroneously assigned by mis- 
take, with the new results (A"''") correctly recomputed in 
this study. Note that these calculations correspond to fc = 1 
(standard treatment of neutral hydrogen collisions). (b) 
Similarly to panel (a), A"^" (filled symbols) and the differ- 
ence <5A (= A""™ — A"'"^; open symbols) are plotted against 
the (non-LTE) equivalent width {EWgo94,)- 
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Fig. 8. (a) Comparison of the carbon abundances (Ac) of 
8 stars derived from tliis study (C I 10683-10691 lines of 
RMT 1) with those obtained by reanalyzing Akerman et al.'s 
(2004) EW data of C I 9061-9111 lines (RMT 3), where 
three different conditions were adopted regarding how neu- 
tral hydrogen collisions (Cjj) are taken into account: (i) LTE 
[k = oo), (ii) standard treatment [k = 1), and (iii) reduction 
of Ch by 1/10 {k = 0.1); note that Ac(i) > Ac(ii) > Ac(iii) 
generally holds. (b) Graphical display of how Ac (RMT 
3) — Ac (RMT 1) varies with a different choice in the treat- 
ment of Ch- 
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Fig. 2. Synthetic spectrum fitting in the 10680-10697 A region, comprising three C I and two Si I lines, for 33 stars in the wide 
range of metaUicity based on the 2009 July data (cf. Paper I). The best-fit theoretical spectra are shown by (blue) solid lines, while 
the observed data are plotted by (red) circles. In each panel (from left to right), the spectra are arranged (from top to bottom) 
in the ascending order of [Fe/H] as in table 1. An appropriate vertical offset (0.05, 0.2, 0.25 for the left, middle, and right panel, 
respectively) is applied to each spectrum relative to the adjacent one. The wavelength scale of each spectrum is adjusted to the 
laboratory frame. 
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Fig. 3. Synthetic spectrum fitting in the 10680-10697 A region for 12 very metal-poor stars (along with Vesta) based on the 
2011 August data (cf. Paper II). The same value (0.05) of vertical offset is applied to each spectrum relative to the adjacent one. 
Otherwrise, the same as in figure 2. 
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Table 1. Parameters of the program stars and the results of the abundance analysis. 
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— 0. 


.04 


+0. 


.15 




HD 115444 


4721 


1.74 


2.00 


— 2, 


.71 


5. 


64 


16 


.9 


— 0. 


,Z4 


— 0. 


,19 




HD 122563 


4572 


1.36 


2.90 


— 2, 


.61 


5. 


80 


23 


. / 


— 0. 


,24 


— 0. 


.13 




BD -1-23 3130 


5000 


2.20 


1.40 


2 


fin 


5. 


99 


36 


.6 


— u. 


9S 

,zo 


-j-U. 


.UO 




HD 140283 


5830 


3.67 


1.90 


— 2, 


.55 


6. 


16 


47 


.0 


— 0. 


,25 


+0. 


.17 




HD 108317 


5310 


2.77 


1.90 


— 2, 


.35 


6. 


32 


60 


.9 


— 0. 


,29 


+0. 


.13 




BD -1-26 2606 


5875 


4.10 


0.40 


— 2, 


.30 


6. 


47 


52 


.3 


— 0. 


,16 


+0. 


.23 




HD 19445 


6130 


4.39 


2.10 


— 2, 


.05 


6. 


53 


62 


.3 


— 0. 


,13 


+0, 


.04 




HD 204543 


4672 


1.49 


2.00 





. / z 


6. 


39 


55 


.0 


n 


,zo 


n 


.4o 




G 29-23 


6194 


4.04 


1.50 


— 1 , 


.69 


6. 


92 


119 


.8 


— 0. 


■ Z 1 


+0. 


.U f 




G 18-39 


6093 


4.19 


1.50 


— 1 , 


.46 


7. 


14 


131 


.9 


— 0. 


,23 


+0. 


.06 




HD 219617NE 


5825 


4.30 


1.40 


— 1 , 


Art 


7. 


13 


104 


.5 


— 0. 


,13 


— 0. 


.01 


North— East component 


HD 219617SW 


5825 


4.30 


1.40 


— 1 , 


.40 


7. 


10 


101 


.2 


— 0. 


,13 


— 0. 


.04 


South— West component 


HD 121135 


4934 


1.91 


1.60 


— 1 , 


.o ( 


7. 


12 


133 


.3 


— 0. 


A A 
,44 


— 0. 


.05 




HD 166161 


5350 


2.56 


2.25 


-1 


.22 


7. 


42 


206 


.3 


— 0. 


.0 I 


+0. 


.10 




HD 194598 


6020 


4.30 


1.40 


-1 


.15 


7 


35 


146 


Q 


— 0. 


,20 


— 0. 


.04 




HD 193901 


5699 


4.42 


1.20 


-1 


.10 


7. 


,39 


118 


.1 


— 0. 


,11 


— 0. 


.05 




HD 201891 


5900 


4.19 


1.40 


-1 


.10 


7. 


,53 


167, 


.1 


— 0. 


,23 


_i_o, 


.09 




HD 148816 


5860 


4.07 


1.60 


-1 


.00 


7. 


,81 


215, 


.3 


— u. 




-j-u. 


97 
.z / 




HD 224930 


5275 


4.10 


1.05 


-0, 


.90 


8. 


,02 


153, 


.7 


Q 


1 1 


_i_n 


38 




HD 6582 


5331 


4.54 


0.73 


-0, 


.81 


7. 


,88 


122, 


.6 


—0 


06 




15 




Arcturus 


4281 


1.72 


1.49 


-0, 


.55 


8. 


15 


73 


.0 


Q 


23 


-j_n 

"t-U. 


16 




HD 142373 


5776 


3.83 


1.26 


-0, 


.51 


7. 


,99 


228, 


.3 


— 


30 


— 


04 




HD 165908 


6183 


4.35 


1.24 


-0, 


.46 


7. 


,96 


243, 


.3 


— u. 


,zo 


— u. 


12 




HD 10700 


5420 


4.68 


0.66 


-0, 


.43 


8. 


12 


148 


.3 


n 


,UU 


1 n 


.Ul 




HD 15096 


5375 


4.30 


0.80 


-0, 


.20 


8. 


70 


246, 


.0 


g 


1 1 


-j-U. 


36 




HD 131156 


5527 


4.60 


1.10 


-0, 


.13 


8. 


,41 


208 


.0 


Q 


09 


i-u. 


00 




HD 141004 


5877 


4.11 


1.17 


-0, 


.01 


8. 


,49 


305, 


.0 


g 


22 


g 


04 




HD 196755 


5750 


3.83 


1.23 


0.09 


8. 


,45 


282, 


.5 


— u. 


,ZU 


— u. 


1 s 

i-O 




HD 161797 


5580 


3.99 


1.11 


0.29 


8. 


79 


302, 


.7 


n 


1 Q 


n 


ri/i 

.U4 




HD 182572 


5566 


4.11 


1.07 


0.33 


8. 


,89 


315, 


.2 


—0 


17 


-f-u. 


02 
















(2011 Au 


gust data) 












CS 30323-048 


6338 


4.32 


1.50 


-3, 


.21 


5. 


,52 


17, 


.6 


-0, 


,28 


+0, 


,19 




HD 126587 


4700 


1.05 


1.70 


-3, 


.16 


5. 


,67 


34, 


.5 


-0, 


,39 


-1-0, 


,29 




G 206-34 


5825 


3.99 


1.50 


-3, 


.12 


5. 


,93 


27, 


.3 


-0, 


,25 


-1-0, 


,51 




G 64-37 


6432 


4.24 


1.50 


-3, 


.08 


5. 


,57 


22, 


.0 


-0, 


,28 


+0, 


,11 


also observed in 2009 July 


HE 1523-0901 


4630 


1.00 


2.60 


-2, 


.95 


5. 


,07 


8, 


.9 


-0, 


,28 


-0, 


,52 




BD -16 251 


4825 


1.50 


1.80 


-2, 


.91 


5. 


78 


35, 


.6 


-0, 


,35 


+0, 


,15 




HD 195636 


5370 


2.40 


1.50 


-2, 


.77 


5. 


,30 


16, 


.8 


-0, 


,33 


-0, 


,47 




G 186-26 


6417 


4.42 


1.50 


-2, 


.54 


6. 


11 


44, 


.3 


-0, 


,19 


+0, 


,11 




HD 186478 


4730 


1.50 


1.80 


-2, 


.42 


5. 


,84 


29, 


.7 


-0, 


,27 


-0, 


,28 




HD 6268 


4735 


1.61 


2.10 


-2, 


.30 


5. 


,63 


17, 


.8 


-0, 


,21 


-0, 


,61 




HD 13979 


5075 


1.90 


1.30 


-2, 


.26 


5. 


,34 


14, 


.8 


-0, 


,23 


-0, 


,94 




HD 221170 


4560 


1.37 


1.60 


-2, 


.00 


6. 


10 


32, 


.4 


-0, 


,21 


-0, 


,44 




HD 216143 


4525 


1.77 


1.90 


-1 


.92 


6. 


,46 


37, 


.0 


-0, 


,17 


-0, 


,16 




Vesta 


5780 


4.44 


1.00 


0.00 


8. 


,54 


283, 


.7 


-0, 


.14 


-1-0, 


.00 


substitute for the Sun 



In columns 1 through 5 arc given the star designation, effective temperature, logarithmic surface gravity, microturbulent velocity dispersion, and 
Fe abundance relative to the Sun, which are the same as adopted in Papers I and II. Columns 6—9 present the results of the abundance analysis. 
Aq is the non-LTE logarithmic abundance of C (in the usual normalization of -Ah — 12.00) derived from spectrum-synthesis fitting, EWioegi is the 
equivalent width (in mA) for the C I 10691 line inversely computed from A^, Aigegi is the non-LTE correction (= — -^^qqqi) for the C I 10691 
line, and [C/Fc] {= —8.54— [Fe/H]) is the C-to-Fc logarithmic abundance ratio relative to the Sun. In each of the two groups corresponding to 
different observational epochs, the data are arranged in the ascending order of [Fc/H]. 
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Table 2. Atomic data of important lines relevant for this study. 



Species 


RMT 


Transition 


A (A) 


Xlow (eV) 




Gammar 


Gammas 


Gam maw 


C I 


3 


3s - 3p ■'Ps 


9061.436 


7.483 


-0.335 


(7.43) 


-5.32 


(-7.59) 


C I 


3 


3s ^P;^ - 3p ^Pa 


9078.288 


7.483 


-0.557 


(7.43) 


-5.32 


(-7.59) 


C I 


3 


3s ^Pg - 3p ^Pi 


9088.513 


7.483 


-0.432 


(7.43) 


-5.32 


(-7.59) 


C I 


3 


3s -^P" - 3p ^P5 


9094.830 


7.488 


+0.142 


(7.43) 


-5.32 


(-7.59) 


C I 


3 


3s 3p° - 3p 3P3 


9111.807 


7.488 


-0.335 


(7.43) 


-5.32 


(-7.59) 


C I 


1 


3s ^'P^ - 3p ^Ds 


10683.083 


7.483 


+0.076 


(7.29) 


-5.40 


(-7.62) 


C I 


1 


3s ^P° - 3p ^Da 


10685.343 


7.480 


-0.276 


(7.29) 


-5.40 


(-7.62) 


C I 


1 


3s 3p° - 3p ^Dt 


10691.240 


7.488 


+0.348 


(7.29) 


-5.40 


(-7.61) 


C I 


1 


3s ^Pa - 3p ^Da 


10707.317 


7.483 


-0.401 


(7.29) 


-5.40 


(-7.62) 


C I 


1 


3s 3p° - 3p ^Ds 


10729.530 


7.488 


-0.401 


(7.29) 


-5.40 


(-7.62) 


Si I 


53 


4p ^Da - 4d 


10689.716 


5.954 


-0.190 


(7.29) 


-4.23 


(-7.29) 


Si I 


53 


4p ^Da - 4d ^F? 


10694.252 


5.964 


-0.060 


(7.29) 


-4.23 


(-7.29) 



All data are wore taken from Kurucz and Bell's (1995) compilation as far as available. RMT is the multiplct number given by the Revised Multiplet 
Table (Moore 1959). In the last three columns are given the damping parameters in the c.g.s. unit: 

Gammar is the radiation damping constant (s~^), log7j-ad- Gammas is the Stark damping width per electron density (cm~'^) at 10"^ K, log(7e/7Ve). 
Gammaw is the van der Waals damping width per hydrogen density (cm~^) at 10^ K, log(7w/ArH). Note that the values in parentheses are the 
default damping parameters computed within the Kurucz's WIDTH program (cf. Leusin, Topil'skaya 1987), because of being unavailable in Kurucz 
and Bell (1995). The meanings of other columns are self-explanatory. 



Carbon Abundances of Metal-Poor Stars from C I 1.068-1.069 /ini Lines 



Table 3. Comparison of the [C/H] results and atmospheric parameters with the literature values. 



Star 








[Ic/HJ 


[G/HJ 


Rcf . 


Line 




(K) 


(cm s } 


(km s ) 


(dcx) 


(dcx) 






BD +44 493 


5510 


3.70 


1.30 


- 3. 68 


-2.85 


{thi.s study) 


C I 10683/10685/10691 




5510 


3.70 


1.30 


-3.68 


-2.37 


ITO09 


CH (4308-4316) 


HD 126587 


4700 


1.05 


1.70 


-3.16 


-2.87 


[this study) 


C I 10683/10685/10691 




4960 


2.10 


1.80 


-2.90 


-2.59 


HON04 


GH (4323) 




4910 


1.85 


2.03 


-2.84 


-2.71 


MGW95 


GH (~4300) 


G 64-37 


6432 


4.24 


1.50 


-3.08 


-2.90 


[this study) 


C I 10683/10685/10691 [2009 data\ 




6432 


4.24 


1.50 


-3.08 


-2.97 


[this study) 


C I 10683/10685/10691 [2011 data] 




6432 


4.24 


1.50 


-3.08 


-3.21 


FAB09 


C I 9061/9062/9078/9088/9094/9111 




6318 


4.16 


1.50 


-3.12 


-2.83 


AKE04 


G I 9061/9078/9094/9111 


BD -18 5550 


^750 


1.40 


1.80 


-3.06 


-2.91 


[this study) 


C I 10683/10685/10691 




4750 


1.40 


1.80 


-3.06 


-3.08 


CAY04 


GH (4224) 




4750 


1.40 


1.80 


-3.06 


-3.08 


SPI05 


GH (4224) 




4790 


1.15 


2.14 


-2.91 


-2.91 


MGW95 


GH (~4300) 




4683 


1.70 


1.70 


-2.87 


-3.26 


MEL02 


GH (4310) 




4750 


1.40 


1.80 


-3.06 


-3.08 


SPI06 


GH (4230) 




4806 


1.72 


1.91 


-2.89 


-3.09 


JOH07 


GH (~4300) 


HD 115444 


4721 


1.74 


2.00 


-2.71 


-2.90 


[this study) 


C I 10683/10685/10691 




4720 


1.50 


1.70 


-2.85 


-3.26 


HON04 


GH (4323) 




4775 


1.68 


2.00 


-2.86 


-3.01 


LAI07 


GH (~4300) 




4721 


1.74 


2.00 


-2.71 


-3.45 


SIM04 


GH (4315) 


HD 122563 


4572 


1.36 


2.90 


-2.61 


-2.74 


[this study) 


C / 10683/10685/10691 




4600 


1.10 


2.00 


-2.85 


-3.22 


CAY04 


GH (4224) 




4570 


1.10 


2.20 


-2.77 


-3.18 


HON04 


GH (4323) 




4600 


1.10 


2.00 


-2.82 


-3.29 


SPI05 


GH (4224) 




4650 


1.40 


2.60 


-2.68 


-3.13 


NOR97 


GH (4323) 




4600 


1.10 


2.20 


-2.62 


-3.02 


AOK05 


GH (4323) 




4600 


1.10 


2.00 


-2.82 


-3.21 


SPI06 


GH (4230) 




4615 


1.27 


2.10 


-2.47 


-3.42 


JOH07 


GH (~4300) 




4600 


1.10 


2.20 


-2.58 


-2.94 


AOK07 


GH (4323) 




4610 


1.32 


2.00 


-2.54 


-2.94 


LAI07 


GH (~4300) 




4572 


1.36 


2.90 


-2.61 


-3.32 


SIM04 


GH (4315) 


BD -1-23 3130 


5000 


2.20 


1.40 


-2.60 


-2.55 


[this study) 


C I 10683/10685/10691 




5224 


2.82 


2.00 


-2.59 


— 2.29 


LAI07 


GH (~4300) 




5285 


2.83 


1.60 


-2.62 


-2.47 


LAI08 


GH (~4300) 


HD 140283 


5830 


3.67 


1.90 


-2.55 


-2.38 


[this study) 


C I 10683/10685/10691 




5849 


3.72 


1.50 


-2.38 


-2.48 


FAB09 


G I 9061/9062/9078/9088/9094/9111 




5630 


3.50 


1.40 


-2.53 


-2.25 


HON04 


GH (4323) 




5750 


3.30 


1.20 


-2.47 


-2.31 


AOK02 


GH (4323) 




5690 


3.69 


1.50 


-2.42 


-2.21 


AKE04 


G I 9061/9078/9094/9111 




5640 


3.28 


1.50 


-2.64 


-2.42 


TOM92 


G I 9078/9088/9094/9111 




5640 


3.28 


1.50 


-2.64 


-2.40 


TOM92 


GH (~4300) 


G 186-26 


6417 


4.42 


1.50 


-2.54 


-2.43 


[this study) 


C I 10683/10685/10691 




6417 


4.42 


1.50 


— 2.54 


-2.65 


FAB09 


G I 9061/9062/9078/9088/9094/9111 




6273 


4.25 


1.50 


-2.62 


-2.41 


AKE04 


G I 9061/9078/9094/9111 


HD 186478 


4730 


1.50 


1.80 


-2.42 


-2.70 


[this study) 


C I 10683/10685/10691 




4700 


1.30 


2.00 


-2.60 


-2.92 


CAY04 


GH (4224) 




4720 


1.60 


2.20 


-2.50 


-2.75 


HON04 


GH (4323) 




4700 


1.30 


2.00 


-2.59 


-2.89 


SPI05 


GH (4224) 




4650 


0.95 


2.71 


-2.58 


-2.86 


MGW95 


GH (~4300) 




4700 


1.30 


2.00 


— 2.59 


— 2.81 


SPI06 


GH (4230) 




4831 


1.78 


1.89 


-2.63 


-2.83 


JOH07 


GH (~4300) 




4598 


1.43 


2.00 


— 2.44 


— 2.96 


blM04 


GH (4315) 


HD 108317 


5310 


2.77 


1.90 


-2.35 


-2.22 


[this study) 


C I 10683/10685/10691 




5234 


2.68 


2.00 


-2.28 


-2.23 


SIM04 


GH (4315) 


BD -1-26 2606 


5875 


4.10 


0.40 


-2.30 


-2.07 


[this study) 


C I 10683/10685/10691 




5910 


3.63 


1.50 


— 2.63 


— 2.35 


TOM92 


G I 9078/9088/9094/9111 




5910 


3.63 


1.50 


— 2.63 


— 2.40 


TOM92 


GH (~4300) 


HD 6268 


4735 


1.61 


2.10 


-2.30 


-2.91 


[this study) 


C I 10683/10685/10691 




4600 


1.00 


2.10 


-2.63 


-3.30 


HON04 


GH (4323) 




4670 


0.75 


2.73 


-2.59 


-3.26 


MGW95 


GH (~4300) 




4705 


1.50 


1.90 


-2.35 


-3.48 


MEL02 


GH (4310) 


HD 13979 


5075 


1.90 


1.30 


-2.26 


-3.20 


[this study) 


C I 10683/10685/10691 




4970 


1.15 


2.35 


-2.65 


-3.25 


MGW95 


GH (~4300) 


HD 19445 


6130 


4.39 


2.10 


-2.05 


-2.01 


[this study) 


C I 10683/10685/10691 




6047 


4.51 


0.80 


-1.94 


-1.99 


GRAOO 


GH (4235, 4365) 




5880 


4.40 


1.50 


— 2.15 


-1.76 


TOM92 


G I 9078/9088/9094/9111 




5880 


4.40 


1.50 


— 2.15 


— 2.20 


TOM92 


GH (~4300) 


HD 221170 


4560 


1.37 


1.60 


-2.00 


-2.44 


[this study) 


C I 10683/10685/10691 




4648 


1.57 


2.22 


-2.15 


-2.70 


BAR05 


GH (4310-4313, 4362-4367) 




4460 


0.75 


1.60 


-2.11 


-2.41 


MELOl 


GH (~4300) [literature value adopted] 




4510 


1.00 


1.80 


-2.20 


-2.89 


IVA06 


GH (~4300) 




4410 


1.09 


1.70 


-2.03 


-2.85 


SIM04 


GH (4315) 




4648 


1.57 


2.22 


-2.14 


-2.71 


ZHAll 


GH (4310) 




4600 


1.50 


1.90 


-1.98 


-2.59 


AOK08 


GH (4323) 


HD 216143 


4525 


1.77 


1.90 


-1.92 


-2.08 


[this study) 


C I 10683/10685/10691 




4360 


0.50 


1.80 


-2.22 


-2.87 


MELOl 


GO (1.56Aim, 2.33Aim) 




4450 


0.80 


2.30 


-2.27 


-2.64 


AOK08 


GH (4323) 


HD 204543 


4672 


1.49 


2.00 


-1.72 


-2.15 


[this study) 


C I 10683/10685/10691 
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Table 3. (Continued) 



[Vol 



Star 






vt 


[Fo/H] 


[C/H] 


Rcf. 


Line 




(K) 


/ — 1 \ 
(cm s ) 


/I — 1 \ 
(km s ) 


(dcx) 


(dcx) 








4570 


1.24 


2.00 


-1.98 


— 2.43 


LAI07 


CH (^-4300) 




4672 


1.49 


2.00 


— 1.72 


— 2.42 


SIM04 


CH (4315) 




4600 


1.00 


2.20 


— 1.81 


— 2.44 


AUKOo 


CH (4323) 


G 29-23 


6194 


4.04 


1.50 


-1.69 


-1.62 


{this study) 


C I 10683/10685/10691 




6194 


4.04 


1.50 


-1.69 


-1.91 


FAB09 


C I 9061/9062/9078/9088/9094/9111 




5966 


3.82 


1.50 


— 1.80 


— 1.70 


AKE04 


C I 9061/9078/9094/9111 


G 18-39 


6093 


4.19 


1.50 


-1.46 


-1.40 


{this study) 


C I 10683/10685/10691 




6093 


4.19 


1.50 


— 1.46 


— 1.45 


FAB09 


C I 9061/9062/9078/9088/9094/9111 




5910 


4.09 


1.50 


— 1.52 


— 1.26 


AKE04 


C I 9061/9078/9094/9111 


HD 219617NE 


5825 


4.30 


1.40 


-1.40 


-1.41 


{this study) 


C I 10683/10685/10691 


HD 219617SW 


5825 


4.30 


1.40 


-1.40 


-1.44 


{this study) 


G I 10683/10685/10691 




5974 


3.88 


1.23 


— 1.48 


— 1.49 


GRAOO 


CH (4235, 4365) 


HD 121135 


4934 


1.91 


1.60 


-1.31 


-1.42 


{this study) 


C I 10683/10685/10691 




4934 


1.91 


1.60 


-1.37 


-1.99 


SIM04 


CH (4315) 


HD 166161 


5350 


2.56 


2.25 


-1.22 


-1.12 


{this study) 


C I 10683/10685/10691 




5270 


2.51 


1.55 


-1.08 


-1.58 


GRAOO 


CH (4235, 4365) 




5350 


2.56 


2.25 


-1.22 


-1.33 


SIM04 


CH (4315) 


HD 194598 


6020 


4.30 


1.40 


-1.15 


-1.19 


{this study) 


C I 10683/10685/10691 




6046 


4.31 


0.80 


-1.08 


-1.08 


GRAOO 


CH (4235, 4365) 




5906 


4.25 


1.30 


-1.17 


-1.07 


AKE04 


C I 9061/9078/9094/9111 




6044 


4.19 


1.00 


-1.16 


-1.13 


SIM04 


CH (4315) 




5910 


4.28 


1.50 


-1.32 


-1.12 


TOM92 


C I 9078/9088/9094/9111 




5910 


4.28 


1.50 


— 1.32 


-1.50 


TOM92 


CH (~4300) 


HD 193901 


5699 


4.42 


1.20 


-1.10 


-1.15 


{this study) 


C I 10683/10685/10691 




5823 


4.58 


0.80 


-1.10 


— 1.21 


GRAOO 


CH (4235, 4365) 




5672 


4.38 


1.00 


— 1.12 


-1.14 


AKE04 


C I 9061/9078/9094/9111 




5750 


4.46 


1.50 


-1.16 


-1.28 


SIM04 


CH (4315) 




5810 


4.83 


1.50 


— 1.22 


— 1.26 


TOM92 


C I 9078/9088/9094/9111 




5810 


4.83 


1.50 


— 1.22 


— 1.90 


TOM92 


CH (~4300) 


HD 201891 


5900 


4.19 


1.40 


— 1.10 


— 1.01 


{this study) 


C I 10683/10685/10691 




5991 


4.30 


1.10 


— 1.13 


— 0.92 


GRAOO 


CH (4235, 4365) 




5909 


4.19 


1.00 


— 1.10 


— 1.04 


aiiViU4 


Cxi (4315; 


HD 148816 


5860 


4.01 


1.60 


— 1.00 


— 0.13 


{this study) 


C I 10683/10685/10691 




5823 


4.14 


1.20 


— 0.73 


— 0.53 


AKE04 


i 90bl/907o/9094/9111 


HD 6582 


5331 


4.54 


0.13 


— 0.81 


— 0.66 


{this study) 


C I 10683/10685/10691 




5408 


4.40 


0.80 


— 0.83 


— 0.79 


iViil/ijUi 




Arcturus 


4281 


1.72 


1.49 


— 0.55 


— 0.39 


{this study) 


C / 10683/10685/10691 




4350 


1.80 


1.70 


— 0.51 


—0.67 


KJ A82 


(^oboUJ 




4286 


1.66 


1.74 


— 0.52 


— 0.09 


KAiVl i 1 


G i oooU/8ooo/9078/91il 


HU 142o7o 


5116 


3. 83 


1.26 


— 0.51 


— 0. 55 


{this study) 


1 10683/10685/10691 




5776 


3.83 


1.26 


— 0.51 


— 0.38 


1 AKOo 


T cnro /coon 

G i ouoz/oooO 


xiD 165908 


6183 


4-35 


1.24 


— O.46 


— 0. 58 


{this study) 


1 10683/10685/10691 




6183 


4.35 


1.24 


— 0.46 


— 0.24 


1 AKOo 


T cnco /coon 

G i ouoz/oooO 


WT^ 1 n?nn 
xiu iUlUU 


5420 


^ . 68 


0. 66 


— 0.43 


— 0.42 


{this study) 


1 lUOoo / luOOD / lUOif 1 




5420 


4.68 


0.66 


-0.43 


-0.20 


TAK05 


C I 5052/5380 


HD 131156 


5527 


4. 60 


1.10 


-0.13 


-0.13 


i^this study) 


C I 10683/10685/10691 




5527 


4.60 


1.10 


-0.13 


-0.14 


TAK05 


C I 5052/5380 


HD 141004 


5811 


4.11 


1.11 


-0.01 


-0.05 


{this study) 


C I 10683/10685/10691 




5877 


4.11 


1.17 


-0.01 


0.02 


TAK05 


C I 5052/5380 


HD 196755 


5150 


3.83 


1.23 


0.09 


-0.09 


{this study) 


C I 10683/10685/10691 




5750 


3.83 


1.23 


0.09 


-0.05 


TAK05 


C I 5052/5380 


HD 161797 


5580 


3.99 


1.11 


0.29 


0.25 


[this study) 


C I 10683/10685/10691 




5580 


3.99 


1.11 


0.29 


0.27 


TAK05 


C I 5052/5380 


HD 182572 


5566 


4.11 


1.01 


0.33 


0.35 


[this study) 


C I 10683/10685/10691 




5566 


4.11 


1.07 


0.33 


0.39 


TAK05 


C I 5052/5380 



For the metal-poor objects with [Fe/H] < -0.8 (i.e., from BD +44 493 through HD 6582), we consulted the SAGA 
database (Suda et al. 2008, 2011) as of the end of 2012 November, in which the [C/H] values are reduced to the 
reference solar carbon abundance of 8.55 (Grevesse et al. 1996). In addition, we also included [C/H] values from 
two papers which were missing in SAGA: 5 stars (HD 19445, HD 140283, HD 193901, HD 194598, BD +26 2606) 
from Tomkin et al. (1992) and 5 stars (G 64-37, HD 140283, G 29-23, G 18-39, G 186-26) from Fabbian et al. (2009; 
results for the 5h = 1 case were adopted and renormalized to the solar abundance of 8.55). Meanwhile, regarding 
stars of [Fc/H] ^ —0.5 (i.e., from Arcturus through HD 182572) where SAGA does not cover, the data were taken 
from Kjaergaard et al. (1982) and Ramirez et al. (2011) for Arcturus, as well as from Takeda and Honda (2005) 
for the other disk stars (an extensive literature survey has not been performed for those non-SAGA stars). The 
reference source is given in the 7th column, and the adopted lines used for C-abundance determination are also 
noted in the 8th column. Our data adopted/derived in this study are expressed in Italic. Stars are arranged in the 
ascending order of [Fc/H] (the metallicity value adopted in this study). 

Key to the references: AKE04 — Akerman et al. (2004); AOK02 — Aoki et al. (2002); AOK05 — Aoki 
et al. (2005); AOK07 — Aoki et al. (2007); AOK08 — Aoki and Honda (2008); BAR05 — Barklem et al. (2005); 
CAY04 — Cayrel et al. (2004); FAB09 — Fabbian ct al. (2009); GRAOO — Gratton et al. (2000); HON04 — 
Honda et al. (2004); ITO09 — Ito et al. (2009); IVA06 — Ivans ct al. (2006); JOH07 — Johnson et al. (2007); 
KJA82 — Kjaergaard et al. (1982); LAI07 — Lai ct al. (2007); LAI08 — Lai et al. (2008); MCW95 — McWilliam 
et al.(1995); MELOl — Melendez, Barbuy, and Spite (2001); MEL02 — Melcndez and Barbuy (2002); NOR97 — 
Norris, Ryan, and Beers (1997); RAMll — Ramirez and AUcnde Pricto (2011); SIM04 — Simmerer et al. (2004); 
SPI05 — Spite et al. (2005); SPI06 — Spite et al. (2006); TAK05 — Takeda and Honda (2005); TOM92 — Tomkin 
et al. (1992); ZHAll — Zhang et al. (2011). 



